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Bifurcation analysis for a free boundary prob- 
lem modeling tumor growth 
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Abstract. In this paper we deal with a free boundary problem modeling 
the growth of nonnecrotic tumors. The tumor is treated as an incom- 
pressible fluid, the tissue elasticity is neglected and no chemical inhibitor 
species are present. We re-express the mathematical model as an oper- 
ator equation and by using a bifurcation argument we prove that there 
exist stationary solutions of the problem which are not radially sym- 
metric. 
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1. Introduction and the main result 



S_h ' Cristini et all. obtained in [5] a new mathematical formulation of an existing 



model (see [3,13,15]) which describes the evolution of nonnecrotic tumors in 
both vascular and avascular regimes. As widely used in the modelling, the 
tumor is treated as an incompressible fluid and tissue elasticity is neglected. 
Cell-to-cell adhesive forces are modeled by surface tension at the tumor-tissue 
interface. The growth of the tumor is governed by a balance between cell- 
mitosis and apoptosis (programed cell-death) . The rate of mitosis depends on 
the concentration of nutrient and no inhibitor chemical species are present. 
This new model is obtained by considering different intrinsic time and length 
scales for the tumor evolution which are integrated by means of algebraic 
manipulations into the model. The model presented in [3, 13, 15], has been 
studied extensively by different authors [1,6-8,13]. It is known that the mov- 
ing boundary problems associated to it are well-posed locally in time [8, 13] 
and, as a further common characteristic, there exists, for parameters in a 
certain range, a unique radially symmetric equilibrium [6-8, 14]. This results 
have been verified to hold true also for the model deduced in [5], cf. [5, 11, 12]. 
The authors of [7-9, 14, 17] show by using the theorem on bifurcation from 
simple eigenvalues due to Crandall and Rabinowitz that in the situations they 
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consider there exists besides the unique radially symmetric equilibrium other 
nontrivial equilibria. Though the problems they consider are different, these 
nontrivial steady-state solutions are asymptotically identical near the circu- 
lar equilibrium. Numerical experiments suggest also for the model [5] that 
there may exist stationary solutions which are no longer radially symmetric. 

In this paper we focus on the general, i.e. non-symmetric situation, 
when the tumor domain is arbitrary and look for nontrivial steady-states of 
the model [5]. Additionally to the well-posedness of the associated moving 
boundary problem, stability properties of the unique radially symmetric so- 
lution are established in [12]. Particularly, it is shown that if G, the rate 
of mitosis relative to the relaxation mechanism is large, then the circular 
equilibrium is unstable, which also suggests existence of nontrivial stationary 
solutions. When studying the set of equilibria we deal with a free boundary 
problem which is reduced to an operator equation between certain subspaces 
of the small Holder spaces over the unit circle h m+l3 (E>). We apply then the 
theorem on bifurcations from simple eigenvalues to this equation and obtain 
infinitely many bifurcation branches consisting only of stationary solutions 
of our model. Near the circular equilibrium these solutions match perfectly 
the ones found in [7-9, 14]. 

The outline of the paper is as follows: we present in the subsequent 
section the mathematical model and the main result, Theorem 2.4. Section 3 
is dedicated to the proof of the Theorem 2.4. 



2. The mathematical model and the main result 



The two-dimensional system associated to the model [5] is described in detail 
in [1 1] . The steady-state solution of the moving boundary problem presented 
there are precisely the solutions of the free boundary problem 
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(2.1) 



The fully nonlinear system (2.1) consists of two decoupled Dirichlet problems, 
one for the rate ip at which nutrient is added to the tumor domain f2, and 
one for the pressure p inside the tumor. These two variables are coupled 
by the fifth equation of (2.1). Hereby «an stands for the curvature of dfl 
and A describes the balance between the rate of mitosis (cell proliferation) 
and apoptosis (naturally cell death). The function / 6 C°°([0, oo)) has the 
following properties 



/(0) = 



and 



/'(VO > for V > 0. 



(2.2) 
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We already know from [11, Theorem 1.1] that 

Theorem 2.1. Given (A, G) G (0, /(l)), there exists a unique radially sym- 
metric solution D(0,Ra) to problem (2.1). The radius Ra of the stationary 
tumor depends only on the parameter A and decreases with respect to this 
variable. 

Hence D(0, Ra) is a solution of (2.1) for all G G M. Thus, we may use 
G as a bifurcation parameter to obtain also other solutions of (2.1). This is 
in accordance with the numerical simulation [5] . 

In order to determine steady-states of (2.1) we introduce a parametrisa- 
tion for the unknown tumor domain Q. Therefore we define the small Holder 
spaces h r (E>), r > 0, as closure of the smooth functions C°°(§) in the Holder 
space C r (§), whereby, § stands for the unit circle and we identify functions 
on S with 27r-periodic functions on R. Furthermore, we fix a G (0, 1) and use 
functions belonging to the open neighbourhood 

V := {p G h i+a (S) : \\p\\ c{ s) < 1/4}, 

of the zero function in /i 4+Q (§) to parametrise domains close to the discus 
D(0,Ra)- Given p G V, we define domain 

fl p := {x G M 2 : \x\ < R A (1 + p (x/\x\))} U {0}, 

with boundary dil p = T p := {R(l + p(x)) x : x G 8} . Given x G T p , the 
real number p{x/\x\) is the ratio of the signed distance from x to the circle 
Ra ■ S and Ra- If il = fi p for some p G V, problem (2.1) re- writes With this 
notation (2.1) is equivalent to the following system of equations 

= /WO m Op, 

Ap = in Q p , 

ip = 1 on T p , 



(2.3) 



p = Kr p - AG— on T p , 



GVi> - Vp - AG| , VN P ) = 



where N p : A(3R A /4, 5R A /A) ffi is the function defined by N p (x) := 
|:r| — Ra — Rap{x /\x\) for all x in the annulus 

A(3R a /4:,5Ra/4:) := {xeR 2 : 3i? A /4 < \x\ < 5R A /4}- 

We re-expressed now problem (2.3) as an abstract operator equation 
on unit circle S. To this scope we introduce for each p G V the Hanzawa 
diffeomorphism Q p : M 2 -> K 2 by 

e p (x)=Rx + ^(\x\~l)p^ 

where the cut-off function ip G C°°(M., [0, 1]) satisfies 

1, H<l/4, 



ip(r) = 



0, |r|>3/4, 
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and additionally max|(//(r)| < 4. In can be easily seen that O p is a dif- 
feomorphism mapping £1 := D(Q, 1) onto f2 p , i.e. O p G Diff i+a (n,n p ) n 
Diff 4+a (R 2 ,IR 2 ). As we did in [12] we define for each p £ V the function 
T(p) := ipoQp, whereby -0 is the solution of the semilinear Dirichlet problem 

| if) = 1 on r p , 
respectively for p £ V and Gglwe set <S(G, p) := p o O p , where p solves 

Ap = in f2 p , 

\x\ 2 ( 2 - 5 ) 
V = kt p -AG^- on T p . 

With this notation, our problem (2.3) reduces to the operator equation 

$(G,p) = in h 1+a {§) (2.6) 

where $:MxV-> h 1+a (S) is the nonlinear and nonlocal operator defined 
by 

$(g, P ) := (gv (r(p) o e; 1 ) - v (s(G, P ) o e; 1 ) - ag|, va p ) o e p . 

(2.7) 

The function $ is smooth $ G G°°(M x V, /t. 1+q (§)), cf. [12], and the steady- 
state D(0, Ra) corresponds to the function p = which is solution of (2.6) 
for all Get 

Therefore, we shall refer to S := {(G,0) : G £ 1} as the set of trivial 
solutions of (2.6). The main result of this paper, Theorem 2.4 states that 
there exist infinitely many local bifurcation branches emerging from E and 
which consist only of solutions of the original problem (2.1). Our analysis is 
based on the fact that we could determine an explicit formula for the partial 
derivative <9 P $(G,0) cf. [12]. Given p £ /i 4 +"(§), we let p = T,kezP( k )x k 
denote its associated Fourier series. Then we have: 



9$(G,0) 



kez, 



J2vk(G)p(k)x\ (2- 



it; 



where the symbol {pk{G))k^i is given by the relation 

_Li fc |3 + _L| W _ G (4!W 



1 1 Au\ u (l) \ 

MG) := --535- 1*| 3 + -dfc| - G --^L- + A /(l) , (2.9) 



and U\u G G°°([0, 1]) is the solution of the initial value problem 

2n+ 1 
r 



-u' = R A f(v )u, 0<r<l, 



u(0) = 1, (2-10) 
u'(0) = 0, 
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when n = \k\ and vq ■= 7~(0). By the weak maximum principle vq is radially 
symmetric, so that we identified in (2.10) vq with its restriction to the segment 
[0, 1]. Particularly, (2.8) and (2.9) imply that 

a(d$(G, 0)) = {^(G) : k G Z}. (2.11) 

Before stating our main result we study first the properties of the symbol 

(fJ-k(Gj)kez- 

Lemma 2.2. There exists M > such that 

* 2TT2 md ^ ± 1 + (2fc + 1)^ + 3) (2 ' 12) 

for all k G N. 

Proof. Let k G N be fixed. Since 

Mr) = 1 + f -Sr f T 2k+1 f(v Q (T))u k (T) dr da, < r < 1, (2.13) 
Jo s Jo 

we deduce that u k is strictly increasing for all k G N. Let now v := u^+i — u/-. 
From (2.10) we obtain 



2fc + 1 , 2 



«' = R 2 f(v )v - -u' k+1 , 0<r<l, 



r r 
w'(0) - 0, 

u(0) = 0. 

Furthermore, we have 

t^o f t^o t AJy y "\2k + 4: 2k + 2 

= -^(^(0)) (2fc + 4) 2 (2fc + 2) < 0, 

which implies by (2.2) that < for t G (0, <5) and some (5 < 1. Thus, 
u is decreasing on (0,(5). Let now t G [0,1] and set m t := minmfi v < 0. 
A maximum principle argument shows that the nonpositive minimum must 
be achieved at t, m t = v(t) which implies u k +i(t) < u k (t) for all t G [0, 1]. 
Particularly, Ufc+i(l) < u k (l). 

We prove now the estimate for u' k (l). Setting M := R 2 a uq{\) max[ 0j i] f'(vo) 
we obtain, due to (2.13), that 

<(!)= f R\ f 'r 2k+1 f'(v (r))u k (r)dr 



r r T 2k+i d T 



Jo 



(2fe + l)(2fc + 3) 



The estimate for u k follows similarly. □ 
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By Lemma 2.2 we know that u' k (l) / u k (l) — > k ^oo 0. Therefore, we may 
define for k 6 N with 

2 u k (l) 

the constant 

Gfc:= A^ ' (2 ' 14) 

2 u fc (l) M j 

which is the only natural number such that fi k (G k ) = 0. We do this since we 
cannot estimate whether [j, k (G), with k small are zero or not. For example, we 
know from [12] that Hx{G) = for all G£K, which makes the things difficult 
when trying to apply bifurcation theorems to (2.6). In virtue of Lemma 2.2 
we have: 

Lemma 2.3. There exists k\ € N with the property that < G k < Gk+i for 
all k>k\. 

Proof. From (2.13) we obtain by partial integration that 

k ■ K(l) - 4 + i(l)) = R\ [ kr 2k+1 f'{v (r))[u k {r) - r 2 u k+1 {r)] dr 

Jo 

T 2k+2 

= R\k ^-~-^f'(vo(r))[u k (r)-T 2 u k+1 (T)} 

1 T 2k+2 



k^-^\f"(v (r))[u k (r)-T 2 u k+1 (r)] 



+ f'(v (r)) K(r) - 2ru k+1 (r) - r 2 u' k+1 (r)} dr 



< i^/'(i)_A_( Ufc (i) - Ufc+xCl)) + £- 



'2k + 2 y Ky ' 1 "(jfe + l)*' 

with a constant L independent of k. Letting now k — ¥ oo we get 

fc-«(l)-4 + i(l)) — ► 0. (2.15) 

k— >oo 

Having this estimates at hand we can prove now the assertion of the lemma. 
For simplicity we set C := 2(f(l)-A)/A > 0,a k := M' fe (l),and6 fe := l+u k (l). 
In virtue of 
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we compute 

Gk+i > G k O 
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oG(3fc 2 + 3fc) > 



a fc+ i(fc-fc 3 ) + a fe (fc 3 + 3fc 2 


+ 2fc) 




(l + 6 fc )(l + 6 fc+1 ) 






a k+ ib k (k - fc 3 ) + a k b k+ i{k 3 


+ 3fc 2 


+ 2fc) 


(I + b k )(l + b k+1 


) 




k 3 (a k - Ofc+i) + a fc+ ifc + a fe 


(3fc 2 + 


2fc) 


(1 + 6 fc )(l + 6fc+i) 






a k+1 b k (k - fc 3 ) + a k b k+1 (k 3 


+ 3fc 2 


+ 2fc) 


(l + 6 fc )(l + 6fc+i) 



Taking into consideration the relations (2.12) and (2.15), we find a positive 
integer k\ such that strict inequality holds in the last relation above for all 
fc > fci. □ 

Finally, we set 



G. 



\k\ - \ kl 
i? 3 1 R\ 



mm 

0<fc<fci 



/(I) _ A 



2 tifc(l) 



A u' k {l) 
2 itfc(l) 



^0 



(2.16) 



The main result of this paper is the following theorem: 
Theorem 2.4. Assume A G (0, f(l)) and that 

2 u (l) 

holds true. Let I > 2 be fixed and k G N, fc > 1 smc/i i/iai Gfc/ > G,. X7ie pair 
(Gfei,0) is a bifurcation point from the trivial solution E. More precisely, in 
a suitable neighbourhood of (Gfc;,0), t/iere exists a smooth branch of solutions 
(G fci (e), p kl (e)) of problem (2.6). Fore — >• 0, we have the following asymptotic 
expressions: 

G kl (e) = G kl +0(e), 

p kl (e)=ecos(kls) + 0(e 2 ). 
Moreover, any other G > G., G (jL {Gfc; : fc > 1}, is not a bifurcation point. 

Condition (2.16) is not to restrictive. Particularly, if Ra = 1 and / = 
ki[o,oo) we have shown in [12] that (2.16) is satisfied. 
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3. Proof of the main result 

The main tool we use when proving Theorem 2.4 is the classical bifurcation 
result on bifurcations from simple eigenvalues due to Crandall and Rabi- 
nowitz: 

Theorem 3.1 (see [2,4]). Let X, Y be real Banach spaces and G(X,u) be a C q 
(q > 3) mapping from a neighbourhood of a point (Ao,Wo) £RxI into Y. 
Let the following assumptions hold: 

(i) G(X ,u )=0, dxG{\ o ,u o )=0, 

(ii) Kcr<9 u G(Ao, uq) is one dimensional, spanned by vq, 

(iii) lm<9 u G(Ao, uq) has codimension 1, 

(iv) d\d\G(\ ,u ) £ Imd u G{X Q ,u ), d\d u G{X , u )v Imd u G{X ai u ). 
Then (Xq,uq) is a bifurcation point of the equation 

G(X,u) = (3.1) 

in the following sense: In a neighbourhood of (Ao,ito) the set of solutions of 
equation (3.1) consists of two C q ~ 2 curves Si and S 2; which intersect only 
at the point (Xo,Uq). Furthermore, Si, S2 can be parameterised as follows: 

51 : (A, u(A)), |A — Ao issmall, u(Xq) = uq, u'(Xq) = 0, 

5 2 : (X(e),u(e)), \e\ issmall, (A(0),it(0)) = (A ,ito), u'(0) = v . 

We want to apply Theorem 3.1 to the particular problem (2.6). As we 
already mentioned $(G, 0) = for all G G M. However, since we can not 
estimate the eigenvalues fj,k(G), 1 < k < fci we have to eliminate them from 
the spectrum of 9$(G, 0), and also to reduce the dimensions of the eigenspace 
corresponding to an eigenvalue (i^ (G) , k > k±, which we may chose to be equal 
to if G is large enough, to one. This is due to the fact that the dimension of 
the eigenspace corresponding to an arbitrary eigenvalue fi^ (G) is larger then 
2, since x and x~ k are eigenvectors of this eigenvalue. 

This may be done by restricting the operator $ to spaces consisting 
only of 2ir/l— periodic and even functions. Given k £ N and / £ N, / > 2, we 
define 

h*+ a (§) := {p £ h k+a (S) : p{x) = p(x) and p{x) = p(e 2 ™ /l x) for all x e §}, 

where for x £ C, x denotes its complex conjugate. Set further V e .i := Vfl 
hf\ a (S). By identifying functions on § with 2n— periodic functions on R we 
can expand p £ /^^"(S) in the following way 

00 

p(s) = J^flfe cos(kls), 

where afc = 2p(kl) for k > 0. With this notation we have: 

Lemma 3.2. Given I > 2, the operator $ maps smoothly R x V Cy i into /ig~J" Q (S), 
i.e. 

0£G°°(]Rx V e j,hl+ a (S)). 
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Proof. The proof is similar to that of [16, Lemma 5.5.2] and therefore we 
omit it. □ 

Finally, we come to the proof of our main result. 



Proof of Theorem 2.4- Fix now I > 2. We infer from relation (2.8) that the 
partial derivative of the smooth mapping $ : Ix V ej ; — > /^"(S) with respect 
to p at (G, 0) is the Fourier multiplier 



d r MG,0) 



OO 

E 

Lfc=o 



afe cos(fcZs) 



Pki (G)a k cos(kls) (3.2) 
is defined by 



for all p = J2^ =0 ak cos(klx) G /i e ~J" a (S), where fj,ki(G), k £ 
(2.9). Our assumption (2.16) implies that p.o(G) ^ 0. 

The proof is based on the following observation: if G > G. and pki (G) = 
then it must hold that G = Gti and kl > k\. Indeed, we notice that if 
/Jfci(G) = 0, for some kl < hi, then 

1_( W )3_ 1 W 

< G„ 



G. < G 



f(l)-A- 



'hi 



(1) 



2 u w (l) 

since, by (2.16) and l>2,kl>2. If kl > h, then p k i(G) = iff 



G 



--L(kl) 3 + \kl 
r a ' R\ 



Au' kl (l) 



G 



kl- 



A-f{l) 



2 u w (l) 

Let 1 < k G N be given such that Gki > G,. From Lemma 2.3 and 
the previous observation we get p m i(Gki) ^ for m ^ k. Our assumption 
(2.16) ensures that the Frechet derivative d p Q(G m i, 0) of the restriction : 
R x V e .i C hf\ a (S) -> ^j Q (§), which is given by relation (3.2), has a one 
dimensional kernel spanned by cos(fcZs). We also observe that its image is 
closed and has codimension equal to one. 

We are left now to prove that the transversality condition (iv 
orem 3.1 holds. Since Gu > G,, our observation implies kl > k\ 
so 

'Au' kl (l) 



of The- 
I- 1 and 



+ A-f{\) \ >0 
2 u k i(l) J 

Further on, we get from relation (3.2) and (2.9) that 

'Au' kl (l) 



d G d p $(0) [cos(kls)] 



+ A - f(l) cos(fcZs) 



2 u w (l) 

and since cos(kls) $ lmd p $(Gki,0) we deduce that the assumption of Theo- 
rem 3.1 are all verified. By applying Theorem 3.1 we obtain the bifurcation 
result stated in Theorem 2.4 and the asymptotic expressions for the bifurca- 
tion branches (G kl {e), p hl (e)) . 
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FIGURE 1. Possible steady-states of the problem (2.1) 

Moreover, if G > G. and G ^ Gki for all k with kl > k% + 1, then it 
must hold that /ifcz(G) ^ for all k £ N, and we may apply [16, Theorem 
4.5.1] to obtain that 9 p <f>(G,0) is an isomorphisms. The Implicit function 
theorem then states that (G, 0) is not a bifurcation point and the proof is 
completed. □ 

The possible steady-states of problem (2.1) are depicted in Figure 1. 
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